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HYBRID  EXPERIMENTAL-NUMERICAL  STRESS  ANALYSIS 


Albert  S.  Kobayashl 
Department  of  Mochanlcal  Engineering 
Seattle*  Washington  98195 


ABSTRACT 

^S>The  hybrid  experl mental -numerical  stress  analysis  technique*  which  saw 
limited  applications  during  the  1950's*  has  been  resurrected  with  the  vastly 
Improved  numerical  techniques  of  the  1970's.  By  Inputlng  the  experimental 
results  as  Initial  and  boundary  conditions*  modem  computer  codes  can  be 
executed  In  Its  generation  and  application  modes  to  yield  results  which  are 
unobtainable  when  only  one  of  the  two  techniques  is  used.  The  hybrid  tech¬ 
nique  thus  exemplifies  the  complementary  role  of  the  experimental  and  numer¬ 
ical  techniques. 


INTRODUCTION 

One  of  the  frustrations  uf  an  experimental  stress  analyst  Is  the  lack  of 
a  universal  experimental  procedure  which  solves  all  problems.  Referred  to  as 
his  second  principle*  Durelll  states  that  "Seldom  does  one  method  give  a 
complete  solution*  with  the  most  efficiency  Cl]".  Examples  of  this  second 
principle  Is  seen  In  photoelastic  coating  and  brittle  coating  techniques  which 
StTO  lonal  strain  gage  testing  In  locations  of  high  stress  concentre- 
f§  aaoaggffcs  Lby  these  two  techniques.  The  hybrid  experimental -numerical 
Ik  4$9l»ss^eJlt!(y^^cechn1que  Is  an  aberration  of  the  above  where  numerical  stKIV*1  *~l 
analysis  remLl  the  second  experimental  method.  -  —  ** 

early  applications  of  the  hybrid  experimental -numerical  stress  «na  l»w1a«aJ  :•** 


analysis  remkp  the  second  experimental  method.  -  —  ** 

early  applications  of  the  hybrid  experimental -numerical  stress 
technique  were  limited  to  separations  of  two-  and  three-dimensional  jpt< 
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stress  In  photoelastic  specimens.  The  first  stress  Invariant  obtained  through 
a  finite  difference  solution  to  the  compatablllty  equation  and  the  maximum 
shear  stress  distribution  provided  by  the  Isochromatics  yielded  the  two  planar 
components  of  the  principal  stresses  [2*3*4],  The  shear-difference  method 
[2,53  and  the  Fllon's  method  [6,7]  used  the  Isochromatlc  and  Isocllnlc  data 
to  Integrate  the  equilibrium  equations  along  a  straight  line  and  a  stress  tra¬ 
jectory*  respectively.  These  single-purpose  numerical  techniques  thus 
provided  only  the  stresses  along  a  specified  Integration  path. 

In  contrast  to  the  above*  the  modem  super  codes  based  on  finite  element 
method*  boundary  element  method  and  finite  difference  method  yield  the  com¬ 
plete  states  of  stress*  strain  and  displacement  for  the  given  constitutive 
relations  and  boundary  and  Initial  conditions.  The  uncertainty  or  the  lack  of 
knowledge  In  these  given  conditions*  however*  limited  the  accuracy  of  the 
otherwise  volumlnuous  outputs  of  these  super  codes.  Inaccurate  numerical 
modeling  procedures  generated  results  with  obvious  errors  and  are  credited  for 
the  resurgence  of  three  dimensional  photoelasticity  In  the  1970’s.  The  hybrid 
experimental -numerical  stress  analysis  technique  of  today  reduces*  If  not 
eliminates*  the  above  uncertainties  In  prescribed  Input  conditions  by  using 
experimentally  determined  boundary  and  Initial  conditions.  The  output  from 
the  otherwise  proven  numerical  techniques  are  either  the  constitutive  relation 
or  the  complete  states  of  displacement*  strain  and  stress  which  cannot  be 
readily  extracted  through  the  use  of  a  single  experimental  technique  In  stress 
analysis.  Thus*  the  hybrid  experl mental -numerical  technique  ts.an  extremely 
efficient  stress  analysis  technique  which  often  provides  more  Information  than 
needed.  The  full  potential  of  the  hybrid  technique*  however*  l»  yet  to  be 
exploited  because  of  the  historic  dichotomy  between  the  theoretlcl ans-turned 
numerical  analysts  and  the  experimentalists. 


*  • 

In  the  following)  the  utility  of  the  hybrid  experimental -numerical  stress 
analysis  technique  will  be  demonstrated  by  some  stress  analysis  problems 
Involving  two-  and  three-dimensional  structural  components*  biomechanics  and 
fracture  mechanics. 


ELASTIC  ANALYSIS  OF  STRUCTURAL  COMPONENTS 

The  numerical  techniques  used  In  modem  hybrid  technique  for  structural 
analysis  are  vastly  superior  to  their  predecessors  since  they  provide  the 
entire  states  of  stress*  strain  and  displacements.  As  a  straight  forward 
extension  of  the  classical  hybrid  technique*  Rao  £81  used  measured  tempera¬ 
ture  and  surface  traction  data  to  solve*  by  the  finite  difference  method*  the 
Beltrami -Michel 1  stress  equations  of  compatlbllty  Interior  to  an  axlsymmetrlc 
solid.  Figure  la  shows  the  end  retaining  ring*  which  Is  shrink-fitted  to  the 
two  ends  of  the  vector  and  which  Is  used  to  contain  the  end  loops  of  rotor 
windings*  In  a  turbo-generator.  The  distributions  of  hoop  stress*  which  Is 
generated  by  shrink  fitting  and  the  centrifugal  force*  obtained  by  the  hybrid 
technique*  three  dimensional  frozen  stress  photoelasticity  and  a  two-dimen¬ 
sional  analog  are  shown  In  Figure  lb.  The  utility  of  this  hybrid  technique  Is 
demonstrated  by  the  author's  quote  of  "The  time-needed  for  the  analysis  Is 
smaller  than  that  required  by  the  time-consuming  and  tedious  shear-difference 
methods"  £8]. 

Figure  2  shows  a  water  turbine  wheel  and  Its  curvilinear  finite  differ¬ 
ence  grid  representation  which  was  analyzed  by  Barlshpolsky  £91.  Frozen 
stress  photoelasticity  was  used  to  determine  the  stress  tensor  on  the  complex 
boundaries.  These  boundary  values  were  Input  to  the  curvilinear  finite  dif¬ 
ference  equations  for  three-dimensional  elasticity  where  the  number  of  equa¬ 
tions  equalled  the  number  of  nodes  and  thus  reduced  the  computational  time  by 
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three  to  six  folds  over  standard  finite  difference  codes.  Die  procedure  Is 
extended  to  steady  state*  three-dimensional  problems  where  measured  surface 
temperature  must  be  Input  In  addition  to  the  measured  surface  tractions  CIO]. 

While  the  specialized  codes  used  in  the  above  hybrid  techniques  are  com¬ 
putationally  more  efficient  by  design*  off-shelf  codes  In  finite  element  and 
boundary  element  methods  are  often  used  for  sheer  expediencies.  For  an 
elasto-statlc  problem*  the  boundary  element  method  Is  more  computationally 
efficient  and  natural  where  the  Input  data  consists  of  experimentally  deter¬ 
mined  boundary  displacements  and  tractions.  When  used  together  with  the 
double  exposure*  laser  speckle  Interferometry*  "the  measured  surface  dis¬ 
placements  become  the  Input  data  needed  In  the  boundary  element  method  to  cal¬ 
culate  the  traction  vectors  at  specified  points  on  the  boundary  til]"  as  well 
as  In  the  Interior  of  the  body.  Moslehy  and  Ranson  till  demonstrated  the 
utility  of  this  hybrid  technique  by  the  excellent  agreements  In  theoretically 
and  experimentally  obtained  stresses  Interior  to  a  cantilever  beam  with  a 
transverse  end  load.  In  a  similar  application  of  the  hybrid  technique*  Bales* 
Sladek  and  Drzlk  [12]  used  the  double-aperture*  laser  speckle  Interferometry 
and  demonstrated  the  advantage  of  the  hybrid  technique  by  analyzing  only  the 
region  of  Interest  of  a  plate-stiffened  frame.  In  this  case*  the  recorded 
displacements  were  input  to  a  simplified  boundary*  which  Is  represented  by  the 
dashed  lines*  of  the  frame  structure  shown  In  Figure  3.  Boundary  elmaent 
method  was  used  to  determine  the  stress  distributions  along  the  three  cross 
sections  shown  In  Figure  4. 

As  a  variation  In  the  above  mentioned  hybrid  technique*  Umeagufewu*  Peters 
and  Ranson  C133  used  two-dimensional  photoelasticity  together  with  a  boundary 
element  code  to  optimize  the  filets  In  a  doubled  notch  tension  plate.  The 
Interior  principal  stresses  obtained  by  the  hybrid  technique  were  used  to 
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more  evenly  distribute  the  load  along  the  not  section  and  thua  resulted  in  a 
batter  understanding  of  the  filet  optimization  prod  lee. 

CORNEO-SCLERAL  ENVELOPE 

The  Interocular  pressure  of  a  human  eye  Is  maintained  at  a  nearly  con~ 
stant  level  of  15  -  20  mmHg  through  a  complex  physiological  system  Involving 
the  mechanical*  biochemical  and  neurological  responses  of  the  eye  C143.  When 
the  outflow  of  the  ocular  fluid  Is  restricted  by  pathological  conditions*  the 
ensuing  Increase  In  Interocular  pressure  eventually  results  In  glaucoma  which 
Is  the  direct  cause  of  13.5  %  of  the  blindness  In  United  States  [153.  Tonome¬ 
try  monitors  this  Interocular  pressure  by  measuring  the  exterior  mechanical 
response  of  the  cornea  which  Is  indented  or  flattened  by  a  tonometer  plunger. 
The  tonometer  reading  Is  thus  affected  by  the  mechanical  response  of  the 
pressurized  corneo-scleral  envelope  which  Is  essentially  a  pressure  vessel 
containing  the  optical  and  neurological  components. 

The  mechanical  properties  of  the  cornea  and  sclera  are  difficult  to  ob¬ 
tain  because  of  the  small  size*  delicacy  and  natural  curvature.  The  commonly 
used  ocular  rigidity  C163#  which  relates  the  pressure  and  volume  of  the 
corneo-scleral  envelope*  Is  a  global  coefficient  and  Is  not  suitable  for 
analyzing  the  local  deformation  process  under  tonometer  loading.  Simple 
tension  testings  of  excised  strips  of  the  cornea  C173  and  the  sclera  C183 
yielded  erroneous  modulus  of  elasticity  and  Poisson's  ratio  by  the  loosening 
of  the  collagen  fibrils  from  the  soft  mucopolysaccharide  at  the  excised  edges. 
In  order  to  overcome  the  deficiencies  of  the  above  global  and  local  approach¬ 
es*  Woo  et  al  [19*203  developed  a  hybrid  experimental -numerical  procedure  for 
determining  the  local  mechanical  property  of  an  Intact  corneo-scleral 
envelope* 
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Woo's  experimental  procedure  consisted  of  Measuring  the  cornea  and  sclera 
deformations  as  well  as  the  volume  changes  of  pressurized  anterior  segments  of 
enucleated  human  eyes.  A  flying  spot  scanner  was  used  to  measure  the  relative 
motions  of  two  white  targets  on  the  cornea  or  sclera  which  were  mounted  on  a 
McEwen-type  chamber  £213.  Woo's  numerical  procedure  consisted  of  matching* 
through  trial  and  error*  the  measured  and  computed  deformations  and  volume 
changes.  A  pressurized  axlsymmetrlc  finite  element  model  of  the  anterior 
segment  of  the  corneo-scleral  envelope  was  used  to  execute  the  finite  elament 
code  In  Its  application  mode  for  this  purpose.  The  resultant  Isotropic* 
trll Inear*  elastic  stress-strain  relations  obtained  for  this  analog  model  of 
the  corneo-scleral  envelope  Is  shown  In  Figure  5.  These  trll Inear  stress 
strain  relations  were  Incorporated  Into  a  finite  element  of  the  total  eye 
which  was  used  to  calculate  the  nonlinear  Intraocular  pressure-volume  rela¬ 
tion.  The  lack  of  bending  rigidity  In  the  cornea  under  the  tonometer  probe 
was  modeled  by  artificially  reducing  the  bending  stiffness  of  the  finite 
elements  In  the  compression  region.  With  this  modification*  excellent 
correlations  between  the  calculated  and  published  experimental  results  were 
obtained  [203. 

The  membrane  shell  elements*  which  were  later  used  to  construct  the 
corneo-scleral  envelope*  shown  In  Figure  6  C223*  removed  the  above  mentioned 
artificial  reduction  In  bending  rigidity  In  the  solid  elements  used  by  Woo. 
Woo's  experimental  data  £193  was  re-evaluated  by  this  membrane  shell  model 
which  yielded  slightly  different  distribution  of  elastic  moduli  along  the 
corneo-scleral  shell.  Such  differences  demonstrates  the  Inevitable  Interde¬ 
pendence  of  the  experimental  data  and  numerical  modeling  of  the  hybrid 
experimental -numerical  technique  where  the  finite  element  model  Is  used  as  an 
analog  model  of  the  experiment  £233. 
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The  anterior  portion  of  the  mwbraoe  flnlto  element  model  was  than  usod 
to  model  tho  dofonaatlon  procoss  undor  tonometer  loading.  Flguro  7  show*  the 
computed  and  measured  [243  relations  of  probe  force  versus  probe  area  under 
applanation  tonometry. 


ELASTIC-PLASTIC  FRACTURE  MECHANICS 

Fracture  parameters  governing  elastic*  elastic-plastic  and  dynamic  frac¬ 
ture*  with  the  exception  of  geometric  quantities  such  as  crack  opening  dis¬ 
placements  and  crack  tip  opening  angles*  cannot  be  measured  directly.  In 
practice*  even  the  above  geometric  parameters  are  difficult  to  quantify  and 
are  often  computed  by  using  analog  models  of  the  crack.  Strain  energy  release 
rate  and  stress  Intensity  factor  In  linear  elastic  fracture  mechanics*  which 
is  a  well  established  analog  model  of  the  crack*  can  be  computed  accurately  by 
using  modern  numerical  codes.  The  various  fracture  packages  for  these  codes 

have  beep  verified  by  a  recent  benchmark  problem  [253  and  thus  should  provide 

/ 

correct  numerical  solutions  to  well-defined  boundary  value  problems.  Once  the 
str4ln  energy  release  rate  or  stress  Intensity  factor  Is  determined*  the  onset 
/of  brittle  fracture  can  be  predicted  If  the  critical  values  of  these  quanti¬ 
ties  are  known;  Their  elastic-plastic  extension*  the  J-Integral*  has  also 
been  used  with  some,  success  In  predicting  the  onset  of  ductile  fracture.  Laws 
governing  other  fracture  phenomena*  such  as  stable  crack  growth  under  large 
scale  yielding*  are  being  Investigated  through  empirical  correlations  of 
fractur*  data  with  computed  fracture  parameters. 

An  approach  which  has  been  used  recently  to  establish  a  stable  crack 
growth  criterion  Is  to  Input  actual  crack  growth  data  as  additional  boundary 
values  to  an  elastic-plastic  finite  element  code.  Kannlnen  et  al.  [263  used 


the  finite  element  code  In  Its  "generation  mode"  to  study  stable  crack  growth 


and  Instability  of  A533-B  steel  and  2219-T87  aluminum  center-crack  and  compact 
spec leans.  A  slellar  approach  was  used  by  Shlh  et  al.  [27]  who  studied  stable 
crack  growth  and  Instability  of  A533-B  coepact  speclaens.  Experimentally 
determined  load-line  displacement  versus  crack  length  relation*  as  shown  In 
Figure  8*  was  used  to  simulate  crack  extension  In  the  two-dimension  finite 
element  model  shown  In  Figure  9.  Two  sets  of  elastic-plastic  analyses  based 
on  J 2  deformation  and  flow  theories  of  plasticity  were  conducted.  Figure 
10  shows  excellent  agreements  between  the  measured  and  computed  applied  load 
versus  load-line  displacement  relations  obtained  by  these  two  numerical  analy¬ 
ses.  The  computed  fracture  parameters  Included  the  crack  opening  displacement 
(COD)*  the  crack  opening  angle  (C0A)»  the  J-Integral  and  the  rate  of  change  of 
J-Integral*  dJ/da.  Since  the  fracture  criterion  for  stable  crack  growth  must 
be  Independent  of  specimen  geometry  and  crack  extension*  these  fracture  para¬ 
meters  were  then  scrutinized  for  constancy  during  crack  extension.  Typical 
dJ/da  and  COA  variations  with  crack  extensions  obtained  by  Shlh  et  al.  are 
shown  In  Figure  11  and  12*  respectively.  Both  Kanlnnen  and  Shlh  concluded 
from  their  hybrid  experimental -numerical  Investigations  that  the  COA  was  an 
attractive  fracture  criterion  for  stable  crack  growth  In  the  presence  of  large 
scale  yielding. 

The  above  studies  demonstrate  the  utility  of  the  hybrid  experimental - 
numerical  technique  In  extracting  candidate  fracture  parameters  which  cannot 
be  obtained  directly  from  either  the  experimental  or  the  numerical  analysis 
alone.  The  hybrid  experimental-numerical  technique  provided  computed  fracture 
parameters*  such  as  J  and  dJ/da*  under  actual  test  conditions  and  not  under 
assumed  test  conditions  which  normally  would  have  been  prescribed  In  pure 
numerical  analysis.  The  technique  also  yielded  numerically  consistent  COD  and 
COA  which  In  theory  are  measurable  but  In  practice  are  difficult  to  determine. 


The  elastic-plastic  finite  element  cod—  with  fracture  paste gm,  Including 
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those  mentioned  above*  are  yit  to  k»  subjected  to  the  rigorous  scrutiny 
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leveled  on  the  elastic  codes.  The  wide  variations  In  the  J-Integrals  of  the 

eld- 1970's  [283  hopefully  have  been  reduced  If  not  eliminated  In  the  elasto- 
plastlc  finite  eleeent  codes  of  today. 


DYNAMIC  FRACTURE 

The  state  of  science  on  dynanlc  fracture  Mechanics  studied  with  dynastic 
photoelasticity  has  been  presented  by  J.  V.  Dally  In  Ms  1979  Nil  Hast  M. 
Murray  Lecture  C293.  He  noted  that  the  crack  tip  state  of  stress  provided  by 
dynamic  photoelasticity  and  dynastic  caustic*  techniques  have  and  will  continue 
to  enhance  our  understanding  on  the  complex  phenomena  of  dynamic  crack  propa¬ 
gation.  Dynamic  fracture  studies  by  these  techniques,  however*  are  limited  to 
photoelastic  polymers  and  to  plane  stress  problems  when  photoelastic  coatings 
or  caustics  are  used.  The  hybrid  experimental -numerical  technique*  when  used 
with  the  generation  node  of  finite  element  or  finite  difference  method  will 
extract  dynamic  fracture  parameters  In  opaque  materials  as  well  as  In  non¬ 
plane  stress  problems.  These  dynamic  codes  which*  unlike  the  well -studied 
static  codes*  required  verlflcalon  prior  to  Its  used  In  dynamic  fracture 
mechanics.  Fracture  dynanlc  results  generated  by  various  two-dimensional 
elasto-dynanlc  finite  dlf foresee  codes  £30*313  and  finite  element  codes 
£32*333  have  been  compared  with  dynamic  caustic  results  of  fracturing  poly¬ 
meric  specimens  £34*353.  Similar  verification  studies  have  been  conducted 
with  dynamic  photoelasticity  £363. 

V-.  *J!  ■  .....  ,  ...  ;■ 

The  verified  numerical  code  can  also  be  used  to  check  ancillary  remits 


•Added  by  the  author. 
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deduced  fro*  the  original  experimental  results*  such  as  the  variation  In  Input 
work*  which  cannot  be  easily  measured*  during  the  fracture  process.  Numeri¬ 
cal  analysis  also  provides  the  transient  energy  partition  for  the  Input 
boundary  and  Initial  conditions.  Such  energy  partition  can  then  be  used  to 
check  the  hypothesis  used  In  deducing  the  experimentally  determined  energy 
partition.  The  legend  of  Figure  13  shows  an  Internally  notched*  semicircular 
photoelastic  specimen  which  was  loaded  with  end  rotation  and  shear  deformation 
[373.  The  reported  dynamic  fracture  toughness  versus  crack  velocity  relation 
[293  was  used  as  a  dyanmlc  fracture  criterion  to  execute  a  dynamic  finite 
element  code  in  Its  application  mode  which  yielded  the  crack  propagation  and 
dynamic  stress  Intensity  factor  histories  [383  which  are  In  good  agreement 
with  the  numerical  results.  Having  verified  the  numerical  modeling  of  the 
photoelastic  experiment*  the  energies  during  crack  propagation  were  computed 
and  plotted  as  shown  In  Figure  14.  The  Internal  consistency  In  the  computed 
energy  partition  verifies  the  basic  postulates  of  negligible  viscoelastic 
damping  and  negligible  energy  dissipation  at  the  finite  specimen  boundaries 
during  the  dynamic  crack  propagation  period. 

A  relatively  simple  application  of  the  hybrid  technique  Is  the  determina¬ 
tion  of  the  dynamic  stress  Intensity  factor  In  an  Impacted  notch  bend  speci¬ 
men.  Measured  time  variations  In  the  striker  load  were  Input  to  the  finite 
element  model  of  a  dynamic  finite  element  code  which  was  then  used  to  compute 
the  time  variations  In  the  dynamic  stress  Intensity  factor  [403.  The  numeri¬ 
cal  code  was  also  verified  by  comparing  the  computed  and  measured  dynamic 
strains  near  the  crack  tip  as  shown  In  Figure  15.  Figure  18  shows  the  vari¬ 
ations  In  dynamic  and  the  corresponding  static  stress  Intensity  factors  with 
time  prior  to  the  crack  propagation.  These  results  show  the  Inadequacy  of  the 
static  stress  Intensity  factor  which  was  computed  by  using  a  static  formula 


and  the  Instantaneous  striker  load.  It  also  Indicates  the  futility  In  Inter¬ 
preting  such  Impact  fracture  response  without  the  use  of  proper  dynamic 
analysis  [41].  As  a  verification  of  codes*  Figure  17  shows  the  agreement 
between  three  Independent  dynamic  fracture  analyses  of  another  Impacted  three 
point  bend  specimen  [42]. 

Figure  18  shows  a  wedge  loaded*  modi  fled -tapered  double  cantilever  beam 

(WL-MTDC8)  specimen  which  was  fabricated  from  plate  glass.  The  specimen  was 

25%  side-grooved  to  guide  the  propagating  crack.  The  flexible*  long  tapered 

beam  sections  was  designed  to  lessen  the  friction  with  the  silicon  carbide 

loading  pin.  The  specimen  was  wedge-loaded  to  fracture  In  a  500-kg  Instron 

testing  machine  and  the  crack  extension  history  was  recorded  by  a  KRAK-GAGE 

and  associated  Instrumentation  [43].  Figure  19  shows  typical  crack  length 

versus  time  data  which  Is  characterized  by  the  unambiguous  Initial  period  of 

crack  acceleration  and  which  has  not  been  observed  In  dynamic  fracture  of 

metals  and  photoelastic  polymers.  The  average  of  the  two  data  sets*  which  Is 

represented  by  a  solid  curve  In  Figure  19*  was  used  to  drive  a  dynamic  finite 

element  code  In  Its  generation  mode.  The  resultant  K^yn  as  well  as  the 

stat 

static  stress  Intensity  factor*  K*  ,  which  was  also  computed  by  finite 
element  analysis*  are  shown  In  Figure  20.  Although  It  Is  not  obvious  from 
Figure  20*  unlike  the  dynamic  fracture  of  metals  and  polymers*  the  crack  never 
arrested  In  these  and  other  ceramics  WL-MTDC8  specimens  [40].  Thus  the  K^yn 
versus  a  curve  In  Figure  21  should  continue  past  the  nominal  static  fracture 
toughness  »  0.73  Mpa  m  as  Indicated  by  the  dashed  lines.  Notable  Is  the 
lack  of  the  typical  gamma-shaped  K<jyn  versus  a  commonly  observe  In  metals  and 
polymers. 


CONCLUSION 


The  hybrid  experimental -numerical  technique  yields  reliable  Information 


which  cannot  be  obtained  by  the  single  use  of  either  the  experimental  or 
numerical  technique.  The  utility  of  the  hybrid  experimental -numerical 
technique  In  experimental  mechanics  Is  demonstrated  by  case  studies  In  two- 
and  three-dimensional  stress  analysis#  biomechanics  and  fracture  mechanics. 
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Fig.  la  End  Ring  Assembly  on  Rotor 
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Fig.  lb  Hoop  Stress  In  End  Ring  Due  to 
Centrifugal  Forces 


Fig.  2a  The  Model  of  the  Working  Wheel 
of  a  Water  Turbine 
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Fig.  2b  Finite  Difference  Network 
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8  Prescribed  and  Simulated  Load-Line  Displacement 
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Fig.  12  Crack -Opening  Angla  vs  Crack  Ext ansi on  for 
4T  Compact  Spaclasn 
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Fig.  13  Stress  Intensity  Factors  and  Crack  Verifies 
of  Internally  Notched  Saul-Circular  Homallte-100 
Specimen  Subjected  to  ^Nota^lon  and  Shear  Formation 
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Fig.  15  Dynamic  Strain  at  Location  (1), 
A533B  Bend  Specimen 


Fig.  16  Stress  Intensity  Factors  of  an  Impacted  A533B  Steel 
Notched  Bend  Specimen  (L  ■  229,  W«51,  B-25,  a«25ram) 
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Fig.  17  Dynamic  Stress  Intensity  Factor  of  a 
Dynamic  Tear  Test  Specimen 
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Fig.  18  Glass  WL-MTDCB  Specimen 
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Fig.  19  Crack  Extension  vs  Time  of  Fracturing 
WL-MTDCB  Glass  Specimen 
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Fig.  20  Stress  Intensity  Factors  and  Crack  Velocities 
in  a  Fracturing  WL-MTDCB  Glass  Specimen 
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Fig.  21  Dynamic  Fracture  Toughness  vs  Crack 
Velocity  Relation  for  Glass 


29 


i 


474 1 Mi  714 1  lab 
?4»474-419 


fart  t  -  9m  rant 

M^mnwlw  «X  mi-«W  Activities 

MIU*  ef  Mval  iMMwk 
UfartMt  af  tha  **7 
Ar  11m  tea,  M  23117 
Act* i  0*4*  474  (2) 

471 

200  _ 


Mrscter 

Office  ef  h*al  hMtrtk 
buck  Office 

444  - - Street 

■setae,  MA  02210 

Of  rector 

Office  of  Hivel  ieeearch 
Iritch  Office 
524  South  Clark  Street 
Chicago,  a  40405 

M rector 

Office  of  Mvel  Baaaareh 
trench  Office 
1050  teat  Green  Street 
Pasadena.  a  91104 

■aval  Ieeearch  Laboratory  (4) 

Code  2427 

Waahlagtoa,  D.C.  20575 

Oefenee  Oocunratatioa  Canter  (12) 
Casaroo  Station 
Alexandria,  Virginia  22514 

Ondereee  Mploaiea  Aeeeercb  Division 
■aval  Ship  taaoarcb  and  Devalopssat 
Ceater 

Her folk  fcval  Shipyard 
Port snout h,  94  23709 
Attn:  Dr.  B.  Falser,  Coda  177 


Navy  (Coa't.) 

■aval  taaearcb  Laboratory 
Week! ag tea,  P.C.  20375 
Attn:  Coda  1400 


1440 

4300 

4390 

4340 

David  V.  Taylor  Bevel  Map  Beeaerck 
aad  Developseat  Canter 

Aaaapolla.  «  21402 
Ac  tat  Code  2740 
24 

241 

■aval  Uaapoaa  Center 
China  Uka.  CA  93355 
Attet  Coda  4042 
4520 

Consaadlag  Officer 

■aval  Civil  Bagiaeeriag  Laboratory 

Coda  L31 

Fort  Be  ■■—a.  a  93041 

■aval  Bar  face  Uaapoaa  Ceater 
White  Oak 

Silver  Spring,  *>  20910 
Atta :  Coda  1*10 
0-402 
Ml 

Technical  Director 
■aval  Ocean  Syita*  Caster 
Sea  Diego ,  CA  92152 

Navy  Oedema  ter  So  aad 
inference  Dlvlalo* 

■aval  Baaaorcb  laboratory 

t.0.  Bos  B537 
Orlando,  FI  32404 


47*i»i714tiab 


■avy  (Coa't.) 

Chief  of  Havel  Operation# 
Oepertseat  of  the  Mevy 
Waehlagtoa,  4.C.  H»S0 
Altai  Coda  04-094 


ftg  (Coa't.) 


d  Dirac tor 
■avid  W.  Taylor  Beval  Ship 
arch  aad  Davalapsaat  I 


Strata*!*  By* tarn 
Dapartsaat  of  chi 
Waahlagtoa,  D.C. 
Altai  tt P-100 


Ftojaet  Offloa 


■aval  AU  Syatasa 
Dapartsaat  ef  the  levy 
Waehlagtoa,  D.C.  20341 
Atta:  Coda  5302  (Aaroepae*  aad  Stroctwree) 
404  (Technical  Ubrary) 

3204  ( Struct area) 

■aval  Air  Develops* at  Canter 
Wars teat ar,  FA  14974 
At  eat  Aerospace  Mechanic* 


Altai  Cede  042 
17 

172 

173 

174 
1400 
1444 

012.2 

1900 

1901 
1945 
1940 
1942 


■aval  Ubdermter  Syatasa  Caster 
■report.  BZ  02440 
4ttai  Dr.  B.  Trainer 


0.1.  ■aval  Aeadasy 
laglaaarlag  Dapartsaat 
Asaapolia,  MD  21402 

■aval  Fee  11  It  las  lag  I  near  I  eg  tf  ad 

200  Stovall  Street 
Alexandria,  VA  22332 
Atta.  i  Cede  03  (Baaaareh  4  Nail  uniat) 
04S 
043 

14114  (Technical  Ubrary) 

Naval  Saa  Syatasa  CoaaM 
Dapartsaat  ef  tbe  Bevy 
Ws eking ton,  D.C.  20342 
Altai  Coda  Oil 
312 

322 

323 
05ft 
321 


■aval  Surface  Weapcae  Canter 
Bahlgraa  Laboratory 
Dahlgraa,  F4  22444 
Acts i  Cede  C04 
620 

Tec be leal  Director 
■ere  Ulead  Naval  Shipyard 
Vallejo,  CA  94592 

O. 4.  Havel  Foatgradnata  School 
Library 

Cade  0S44 

Ntataray,  CA  93940 

Mibb  Institute  ef  Naval  Arehltacter* 
At tat  Librarian 
Crane eat  Beach  load,  OLaa  Cove 
LOU  Ulead.  WT  11542 

4ns 

Oassaadiag  Officer  (2) 

D.S.  Amy  Be  search  Office 

P. 0.  Bra  12211 

Beeaerck  Trte^le  Fork,  NC  27709 
Altai  m.  J.  J.  Mar  ray,  CSD-AA-IP 


474 1 HP: 716  slab 
74nd  74-dlP 


474tMFs 714i lab 
74v474-#l9 


Away  (Coa't.) 


Air  Force  (Coa't.) 


Watervliet  Araeaal 
NAG0S  Be anarch  Ceater 
WatervlUt,  Wf  12149 
Altai  Direetav  of  leaeereh 


Chief  Applied  Maehaaice  Oreap 
O.s.  Air  Force  IBatltut*  ef  Technology 
tfrlfbc-Pacterse*  Air  Urea  Base 
Dayton.  <*  45433 


0.1.  Any  Notarial#  and  Beh 
Baaaareh  Center 
Water tone.  MB  08172 
Attat  Or.  B.  Shan,  ODS-T 


0.1.  Arsy  Missile  Be  search  and 
Develepnsat  One ter 
Badstene  Bclaetlflc  Isfecsatles 


Center 


Chief.  Deems* 
Bed at one  hr  easel 


Amy  Bveeersb  end  Daralepsaat 
Caster 

Pert  Balvelr ,  FA  14040 


Chief,  Civil  Hagiaearlag  brasah 
WtBC,  Ba search  Dlvlalo* 

Air  Fore#  Weapon*  Laboratory 
Klrtlead  Air  Farce  Be a* 

Albueergsa,  W  47117 

Air  Force  Office  ef  Selestiria  Baeearch 
Bell lag  Air  F»rc*  Baeo 
Waehlagtoa,  D.C.  20332 
Atta i  Maehaaice  Division 

Dapartsaat  ef  th*  Air  Fare* 

Air  Oalverslty  Library 
Mama  11  Air  Fares  Base 
Maatgasary,  At  Mill 


Other  Give resent  Activities 


Mtleeal  Aeraaaatiea  end  Sgeae 
Bdstaletretlae 

Bersttaree  Heeaarah  givlalaa 
Lam  ley  — «*■>  Smear 
leagley  Steeles 
■anptea,  MB  SIMS 

Mtlcsel  Aatsaaatiu  and  Spaea 
Adaiaiecrstfcea 

AeaeaUte  Adaiaistrstar  for  Advanced 
Waehiagtea,  ».C.  209*4 


Air  Ferae 

Mright-Fitterece  Air  Faroe  Sect 
laytee,  m  *5433 
Altai  AfFBL  (fl) 

(VW> 

(PM) 

(my 

AIM.  (MM) 


Chief,  Tbatiag  end 
0.4.  Genet  Oeerd 
1)00  B  Street,  BW 
Maebiagtea,  D.C. 

Technical  Mreeter 
Marine  Corps  Develop 
aad  admit lra  Case) 
Qneat  ioe,  VA  281)4 


Director  Hafesm  Hnccrck 

am  Bagleeerlm 
Tae  bales  1  Ubrary 
Bees  3C12B 
Tbe  PmtefBN 
Maebiagtea,  D.C.  20MI 


Bevclepmae  DivUies 


Other  Coverraaat  Actlvltlaa  (Coa't.) 
Dr.  H.  Gaea 

Matleasl  Sc  lanes  fomdatlea 
lavlromeseal  Ba  search  Division 
Waahlagtoa ,  D.C.  20550 

Ubrary  ef  Co  agrees 

•cleaee  and  Technology  Division 

Washing toa,  D.C.  20540 

Director 

Dafanaa  Maclaar  Agaaey 
Washing ton,  D.C.  20)05 
Atta  i  SMI 


FAIT  2  -  Contractors  and  other  Technical 
Collaborator* 

Daivaraitiea 

Dr.  J.  Tinsley  Odea 
University  of  Texas  at  Aeatla 
M5  Bag  invar  lag  Belaace  Bnildlag 
Asst in,  TX  74712 

Frefaasar  Jallna  Miklosita 
California  last l tat*  ef  Technology 
Division  of  Mag laser lag 
and  Applied  Sciences 
Pasadena,  CA  91109 


Mr.  Jarose  Par  ah 
Staff  Specie! let  far  Mteriele 
ead  It  me  tarn* 

0044BAB,  The  Faatagoa 
Haas  30X049 

Waahlagtna,  D.C.  40MI 

Chief,  Airfrma  ead  tgefpsmf  Branch 

Ft-120 

Office  ef  Plight  ftendaide 
Federal  Aristlas  Agency 
Waahlagtaa,  D.C.  2055) 


Or.  Bare  Id  Uehesite,  Dam 
Ob  he*!  ef  Oigiaeerlag  end 
Applied  Bciesca 
Oaerga  VuklagtM  Oaiveraity 
Maebiagtea,  i.C.  20052 

Prefaeeer  111  Dtambarg 
Callfersie  Xnetltnte  of  Tecbsalagy 
Mvlelea  af  Mglacarlag  and 
Applied  Bclescea 
Faaadaaa,  CA  91109 


Mtiesal  An  ad  any  ef  Ft  lessee 
Mtlesel  Baaaareh  Oaeseil 
Ship  Bell  Baaaareh  Gamuts* 

2101  Orastltstlea  Arens* 

Meahiagte*,  D.C.  H41« 

Altai  Ik.  A.  L  Lytle 

Mtlesel  Bciesea  Faasdatlea 
■mtAurlm  Marhaatea  Baa  ties 
Mviales  af  bglssis 
Meahiagte*,  D.C.  14550 

Fieatiasy  Araeaal 

Flattie*  tbe  laical  Oeelmtlca  Oraaar 
Attai  Taahaieel  Ufematiaa  leet&es 
lever,  HI  07410 

Mar  it  Isa  Amulet  reties 
Office  ef  mrltftsa  Tbctselmp 
Ucb  eM  OmetlCttlm  dm.,  W. 
Maehiagtea,  D.C.  402)0 


University  ef  Call Ural* 

■•garment  ef  Mechealeal  lagtaaerlag 
tavhaWy,  GA  94720 


Frstaeear  A*  J.  Oamili 
Oakland  Oaiveraity 
■eheel  ef  Biglmeflm 
Beekaater,  HD  4404) 


Fnfamar  V.  L.  M Haggle 
Oalemii  Oaivsreiky 


Hepsrmmt  af  CivU 
Man  Verb,  HI  14017 


Tha  MmMty  ef  Uvsrpasl 
Bmanseat  ef  I 
P.O.  bra  147 
Bmmlss  Mil 
Uvsrpeal  149  JK 


47«iVtniiM 


4 


fctWflttM  <Cm'C.) 

InUnw  J.  ttimt 
hlytMlMte  iHtUM  «f  Mm  lNk 
Depart neat  «f  tachanlcal  aad 
tarMiwi  b|iM«rii| 

WJ  Jay  kmt 

•r  oak lye,  m  11201 

>r*l •  I.  4.  IcHaptry 
TV***  AM  lhl*«r<lt|r 
(iMrtMiit  of  Civil  h|iBMrlm 
Col las*  Uatle*.  TX  77443 

Professor  Walter  0.  Pllhsy 
University  of  Virginia 
ha  search  Lahar  stories  for  the 
Baplnevring  Science*  mi 
Aral 1*4  Science* 

Charlottesville,  IK  22401 

Profveaor  K.  0.  W1 1  Inert 
Clarkson  Coll»av  of  tachaolcgy 
teperteviu  of  Muchaalcil  Koptneertei 
Poeedea,  Iff  13474 

Dr.  Wat tor  B.  Hauler 
Teen*  44  M  University 
Aeroepac«  Engineering  Us  pec toast 
Col lag*  5tatloa,  TX  7714) 

Dr.  Hub  win  4.  Kta«l 
University  of  Kr  1  iom 
Depertaent  of  Aerospace  a*t 
Mechanical  Knot ovaries 

Tucaon,  Kf.  4)721 

Dr.  S.  J.  Fvnvee 
Carnegie- Halloa  University 
Oapertavnt  of  Civil  B<*ln**rlag 
Schentey  park 
Pittsburgh.  FA  1)212 

Dr.  Bnnald  l.  Huston 
Depart vent  of  Engineering  Analysis 
Oaivvrslty  of  Cincinnati 
CUciaoati,  QM  4)221 


Fry fa a a* r  «.  c.  V.  Ilk 
UM|I  Be  l vet city 
UetUat*  *f  Praetor*  mi 
1*1  id  Mechanics 
Bathlahaa,  ft  ltDIS 

Pro  fa  **o  r  Alkart  I.  taheyaahl 
Univvralty  of  Washing  toe 
Dspartoant  of  Mechanical  toslaatrii 
Soattla,  UA  4B10) 

Frefasaor  Oaaiat  Frederick 
Virginia  Pol ytoc Hole  1  not  1  tot*  Mi 
State  University 

0a par toe* t  of  Baslo**riat  Nee key ice 
Blacksburg,  VA  UNI 

Professor  A.  C.  tries** 

Prlacato*  University 
Depart** at  of  Aeteepaa*  mi 
Mechanical  Science* 

Prise# ton,  t u  0«)40 

Professor  |.  M.  tee 

Staaferd  Uatvaralty 

Divielea  of  Eag  lea*  ring  Hack** tea 

Staa ford ,  CA  4420) 

Professor  Alkart  f.  Mag 
Wayne  State  Univereity 
•lease kanice  Bvaearek  Center 
tatmlt,  Ml  4 4242 

Dr.  V.  I.  Modaaon 
Usyna  State  University 
School  of  Nedlclne 
Detroit ,  HI  44202 

Dean  •.  A.  Bolvy 
Herthunater*  Unlvsrslty 
Depertaent  of  Civil  Engineering 
tv east on,  IL  40201 


Pv*f*o**r  P.  6.  Jr. 

Bnlvorsity  *f  HI ana seta 
•sparraaat  *f  Sara apace  Daglaearlag 
•ad  Mac has lea 
Wlaaa  spa  I  te,  m  IMS) 

Br .  p.  c.  tracker 
Mai  varsity  of  llllanls 
tea*  of  Baglaaarlag 
Urbane,  IL  41*01 

Prefvssor  N.  J.  Ms v ark 
University  of  til  las la 
Departaeat  of  Civil  Baglaeerlag 
tfrkau,  IL  4ID0) 

Frefasaor  K.  ketsaaar 

Del  varsity  ef  California.  San  Map* 

B*pert*e*t  *f  Applied  Machs* Ice 

U  Jolt*.  CA  42017 

Professor  Hi 11  tan  A.  Mesh 
Univereity  of  Nanaaehnaatte 
Departaeat  ef  Hachaeice  a*4 
Aerospace  Baglaaarlag 
Asher  at,  MA  01002 

Profaaaar  C.  Hermann 

Stan ford  Uni varsity 

Departaeat  of  Applied  Mechanic* 

Stanford .  CA  44)0) 

Profaaaar  J.  0.  Ac ha* tech 
tart hue  at  Ualvaralty 
Departaeat  ef  Civil  BAglaaerlag 
Ivaaatoa,  1L  40201 

Prafoaaar  I.  B.  Dong 
Ualvaralty  of  California 
Departaeat  of  Mechanic* 

Us  Angeles ,  CA  40024 

Professor  Burt  Paul 
Unlvsrslty  of  Paaaaylvsaia 
Teona  School  of  Civil  sad 
Hachaateal  Eaglaaarly 
Philadalphie.  PA  14104 


tT*t»i)une 

IWM-UI 

UwnlitM  (Cm'i.) 

fntHw  i.  ».  u. 

irrKMt  Wnnii) 

Bapart— H  of  Chaalcat  Cngtaaertag 

mi  Hmuaaif 

Syracuse,  m  12210 

Pmfnsoar  |.  Bodnar 
Tech* lea  Ml  Paundetfen 
Half*.  Ur  eel 

Profnssar  Wnrnvr  Caldaalih 
University  of  California 
Depart sent  of  Nucheafrel  Cnplnnur <M| 
Berkeley.  CA  44720 

Professor  a.  &.  klvlln 
Lehigh  University 
Center  for  Application 
of  Netheaetlca 
Both! shea,  PA  1401) 

Professor  f.  a.  Coasaralll 
State  University  of  Boa  York  at 
Buffalo 

Division  of  Interdisciplinary  Undine 
Here  Parker  Re*  1  oeer  la*  Bui  Min* 
Chanietry  toed 
Buffalo.  MV  t  "14 

Professor  Joseph  V.  B»a*> 

Orexel  University 

taper taent  of  Mnchat>«*»l  Engine-Meg 
and  Mechanic* 

Philadalphie.  PA  14104 

Profosaar  B.  C.  tees  Maun 
Ualvaralty  af  Maryland 
Aerospace  Engineer  lag  Ch- par  revet 
Cal lap*  Park,  ND  24742 

Professor  Joseph  4.  Clerk 
Catholic  Univereity  of  Aaertee 
Departaeat  of  Nor heniml  Eng t near  ( a* 
Weehl niton,  D.C.  20044 


Univerettien  (Coa't.) 

Ur.  Saauel  B.  Bat  dor  f 
Ualvaralty  af  California 
School  of  Baptnearian 
end  Applied  Science 
La*  Assets*.  GA  40024 

Professor  Isaac  Pried 
Baato*  Ualvaralty 
Dapnrtnont  of  NathaaaUcs 
•eaten,  m  0021) 

Profaaaar  B.  K ran  pi 
Banaanlanr  Polytechnic  laatitut* 

Vt vista*  af  lap i near ly| 

Bap  laser  lap  Mm  hew  Ira 

Tray,  M  121B1 

ta.  Jack  B.  flnaaa 
Ualvaralty  af  Balauar* 

Departaeat  af  Mat  Man leal  and  Aaraapnm 
Bnptnaorinp  and  tha  CMntpr  far 
Canpaalt*  Hater  lata 
Msunrh,  ■  mil 

Mr.  J.  Daffy 
Oman  University 
Bivlalaa  af  Baitapactgg 
Prevldanae ,  Bt  02412 

Br.  J.  L.  Dead  tea 
Carnap te  Malian  Mil varsity 
taperkaaaf  af  Msahaalaal  Oao  inner  uo 
PtCtsDnrph,  PB  12912 

Dr*  »*  *.  Parade* 

Okie  Btata  MU  varsity  BsasaraB  MUM 
Be  park— at  af  Mala  ear  lap  UmBwIm 


Ualvaralty  of  taaeoylvaet* 
UsypTBumB  of  MBtallamy  nod 
Motor  la  1  *  Miens* 

•all#*  af  U# laser M|  and 
Applied  fcUeat 

fMiutiuu,  pa  lti* 


474 i MP: 714: lak 
7Bu4 74-414 

Universities  (Can't.) 

Br.  Jackson  C.  S.  tent 
Ualvaralty  af  Maryland 
Department  of  Arc  hen I  cel  lay laser  lap 
Collnpa  Perk,  m  20742 

Pvofeaaor  T.  f.  Chant 
Ualvaralty  af  Akron 
tapereavat  of  Civil  Baplaeerina 
Akron,  OH  44 W) 

Profaaaar  Chariot  W.  tart 
University  af  OMahaae 
School  af  Aaraapacs,  tachaaleal, 
aad  Unclear  tap laeerlap 
Mr* an,  OK  7)014 

Dtafsasar  tat  ye  n.  Atluri 
Oasrpte  Inst  Unto  af  Tsahaalapy 
Mhnal  af  Mipioaarlnp  aad 
MthMlt* 

Atlanta,  00  20)22 

ffvpfpppar  Orehna  y.  Oaray 
Onivaraity  af  Tanas  at  Aon* in 
tapartaant  of  Aaraapaaa  B%l  gear  lag 
and  iagiaoefiag  iMchnnlaa 
AaPtM.  tl  70711 

Bp,  g,  |.  Mono 

•alvarslty  af  llllnnt* 

Oapnrtnant  af  taeamtical  and 
Applied  IN* has lee 
BPOaaa,  u  dtdoi 


Bum  AviOyna 
Oinalaa  af  tmm 
•stance*  oar  parse  MA 
■wn  tat  u*.  nb  oi  an 

Aryans  Hit  tonal  laboratory 
Uhrary  Barvloaa  Oaparf  at 


474:WP«7l4:Uh 

70*474-414 


*rch  fast I tu to*  (Coa’t.)  Industry  and  1 


tk  Inatiiataa  (Can't*] 


Or.  N.  C.  Jasper 

Canheldge  Aceuatleel  Associate* 

H  * lady*  Avon**  Du tansies 
Ceabrldge,  HA  02144 

Or.  P.  tad la* 

Oanaral  Oynnalt*  Oar para ties 
Bloc  trie  Boat  Mvlelen 
Crs toe,  CT  04)40 

Or.  J.  I.  Oraonapaa 

J.  6.  lap  laser  tag  Be  anarch  Inane  is  toe 

2B21  Mania  Drive 

Balt tnere,  IS  1121) 

Hi  apart  Mtua  ahlphullplnp  Mi 
Ory  tack  CM  pa  ay 
Uhrary 

■depart  tan.  «A  2M41 

•r.  if.  r.  Baatch 

Me Beans 11  Basil**  Carper a tie* 

5201  Only  Avenue 

Haiti ap tan  Batch,  QA  42047 

Or.  m.  m.  ihranaaa 

Oast has al  Oaaanrah  laptitata 

0)00  Calekra  Bead 

•an  Antoni*.  Tl  70204 

•r.  B.  C.  tatart 
Baathaaat  tataarah  laatttaha 
$m  Calahr*  Bead 
Baa  Antonia,  ft  71204 

Dr.  N*  l.  Daraa 
Paid  linger  Asset  Utes 
110  BaM  94th  It  rent 
He*  Tark.  nr  10012 

Dr.  T.  L.  taare 

laakhaad  Msallaa  and  Bpnee  Oanyay 
M)l  M  navy  it  mat 
Dal*  Alu,  04  44 M4 

Nr.  VI 11  las  Caywnd 
Apt  Had  Thyalr*  Uktratary 


Or.  Oahnrt  K.  taw hen 
Fee if  ice  Technology 
P.0.  Ban  140 
0*1  Mar,  CA  BB014 

Or.  N.  D.  tannines 
Battalia  CaIMkat  Lahore  ter 
V»  Ring  Avenue 
Celunhua,  00  42201 

Dr.  A.  A.  Bachratn 
Daadnlaan  Aaaacletaa,  Inc. 
tprlnglaha  ta search  Band 
1)110  Trader  tch  Band 
Wnsdhta*,,  m  21747 


tanm  Barvlna  Carpet* ties 
r.o.  tan  Mi) 

Bent lag  tan  Bench,  CA  4BB44 

N.  Oahnrt  «.  OUtell 
Applied  Oeiaaoo  and  ttahnalagy 
2244  Berth  tar  my  TIM  Chart 
Wits  DO 

U  Jolla,  CA  «2M9 

ta.  Bavin  Mums 
taatlaghaasa  Roe  trie  Carp* 
Atvaacad  Bane  tar*  Mvinlaa 
D.0.  Bm  19B 
Mad  loan,  UA  1MB) 

ta.  taraaid  Shoffar 
talytashnla  laatlhaaa  af  tav  tar* 
•apt.  of  taahaalaal  had  Anraaphea 
Dag ins a r lap 
)))  Jap  Btmat 
tanpkiyn,,  VF  Uttl 


security  ci Attir icatiom  or  tm»  page  n«<  .»><*», «<o 


REPORT  DOCUMENTATION  PAGE 


UWA/DME/TR-83/47 


«.  TITlt  fan*  Submit) 

Hybrid  Experimental -Numerical  Stress  Analysis 


T.  author^ 

A.  S. 'Kobayashi 


>  PERFORMING  ORGANIZATION  name  ano  aooress 

Department  of  Mechanical  Engineering  FU-10 
University  of  Washington 


it.  controlling  orncc  name  ano  aoorcss 

Office  of  Naval  Research 
Arlington,  VA  22217 


4  MONITORING  AGENCY  NAME  a  AOORCSSftff  dlliatant  from  Contrasting  Ot/lc •) 


14  DISTRIBUTION  STATEMENT  (at  thi a  Report) 


Unlimited 


IT.  distribution  STATEMENT  (•!  ill,  sb.it.el  tnft.d  In  Block  20,  It  dllltrtnl  hum  MtpoH) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


».  ktClkllliT'l  CAT AlOG  NUMBER 


t.  TYPE  OF  REPORT  I  RERIOO  COVEREQ 


Technical  Report 


N0001 4- 7 6- C- 0060 


PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

NR  064-478 


It.  REPORT  OATS 

April  1983 


It.  NUMBER  or  PAGES 

30 


IS.  SECURITY  CLASS.  f»f  (Ala  npM) 

unclassified 


Urn.  declassification/ DOWNGRADING 

schcoule 


If.  KEY  WOROS  (Cwillnua  an  ravaraa  aid*  It  nacaaaary  and  Identity  If  llech  nunlar) 

Finite  element  method.  Boundary  element  method,  finite  difference  method, 
biomechanic,  fracture  mechanics 


10.  ABSTRACT  (Canelnua  an  raaaeao  aldm  It 


•ary  and  Identity  by  black  number) 


The  hybrid  experimental -numerical  stress  analysis  technique,  which  saw  limited 
applications  during  the  1950's,  has  been  resurrected  with  the  vastly  Improved 
numerical  techniques  of  the  1970's.  By  1rtputrng~ttTe~expeHfli€!rM,rfiisuUs 
as  Initial  and  boundary  conditions,  modern  computer  codes  can  be  executed  1m 
Its  generation  and  application  modes  to  yield  results  which  are  unobtainable 
when  only  one  of  the  two  techniques  Is  used.  The  hybrid  technique  thus 
exemplifies  the  complementary  role  of  the  experimental  and  numerical  techniques 


DO  |j*£*T>  1473  EDITION  or  I  MOV  ••  IS  OBSOLETE 

*/R  BIO 2*0 1 4*  4AOI  | 


Unclassified _ ' 

BSCURM  Y  CLAOOlPlCATlON  OP  THIB  PAOE  (Rm  i 


